Aharonov-Bohm oscillations in the local density of topological surface states 
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We study Aharonov-Bohm (AB) osciUations in the local density of states (LDOS) for topological 
insulator (TI) and conventional metal Au(lll) surfaces with spin-orbit interaction, which can be 
probed by spin-polarized scanning tunneling microscopy. We show that the spacial AB oscillatory 
period in the total LDOS is a flux quantum ^o=hc/e (weak localization) in both systems. Re- 
markably, an analogous weak antilocalization with $o/2 periodic spacial AB oscillations in spin 
components of LDOS for TI surface is observed, while it is absent in Au(lll). 
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Spin-orbit interaction (SOI) effects in semiconductors 
and metals have been an active theme in the modern 
condensed matter physics, but especially, the topological 
insulators (TI) , which have been detected in series of two- 
dimensional [iHll and three-dimensional 04111 materials, 
suggest new directions for this field since the extraordi- 
narily strong SOI exists in TI. The helical spin structure 
of electrons in TI acquire a spin-orbit induced nontriv- 
ial Berry phase of tt after a 27r adiabatic rotation along 
the Fermi surface. This Berry phase corrects the quan- 
tum constructive interference between a closed trajec- 
tory that the electron passes and its time-reversal coun- 
terpart, and thus can give rise to the weak antilocaliza- 
tion (WAL) signature in the quantum coherent magneto- 
transport coefficients. Many efforts have been devoted 
to observing WAL in the HgTe quantum wells [H, [11] , 
Bi2Te3 [mill, BizSea and other materials. Very 

recently, for example, in a transport experiment [17] of 
Bi2Se3 nanowires, the Aharonov-Bohm (AB) oscillation 
with anomalous period ^Q—hc/e (twice the conventional 
period) of magnetoconductance was observed, while the 
WAL induced $o/2 period was absent. The period of 
these oscillations in conductance is determined by the 
doping level and the disorder of the TI nanowire [H, [ll] . 
In addition, an energy gap at the Dirac cone opened by 
magnetic doping in TI films could also induce crossover 
from the WAL to the weak localization (WL), which is 
tunable by the Fermi energy and the gap [23 ]. 

In this letter, for exploring the analogous WAL in topo- 
logical surface state, we consider an imaginary AB in- 
terferometer consisting of a spin-polarized scanning tun- 
nehng microscope (SP-STM) tip [25l - [27| and two iden- 
tical nonmagnetic impurities apart tens of nanometers 
on the TI surface as shown in Fig. [TJ When an elec- 
tron travels along the clockwise and anticlockwise loops 



(r^ri^r2^r) enclosing a finite area, the quantum in- 
terference phenomenon occurs. The interference contri- 
bution to the local density of states (LDOS) is affected 
by an external magnetic field B via the AB effect arising 
from the threaded magnetic flux [28]. We will show that, 
on one hand, the AB oscillatory period of the total LDOS 
ANl {r,uj,B) (the subscription L represents the loops en- 
closed by the scattering paths of the surface electrons) is a 
flux quantum $0j which is an analog of WL phenomenon. 
On the other hand, the strong SOI in TI materials can 
modify the electron AB interference effect via the chi- 
ral spin rotations during the non-collinear multiple scat- 
tering processes, and thereby the analogous WAL effect 
with $o/2 period in AB oscillations will be observed if 
the spin-resolved LDOS [AnI {r,uj,B) and AiV^ {r,uj,B)] 
are measured. Furthermore, for comparison, we briefly 
discuss the AB effect in the LDOS of conventional metal 
Au(lll) surface with weak SOI, and find that the anal- 
ogous WAL phenomenon is absent in Au(lll) surface. 
Our finding may provide a useful method to characterize 
the topological surface states. 

We describe the TI surface, on which two nonmagnetic 
impurities are adsorbed, by a low-energy effective Dirac 
Hamiltonian written as 



H = va- ■ {z X q) + V (r) 



(1) 
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with the Fermi velocity w=333 (287) meV-nm in Bi2Se3 
(Bi2Te3) [29j . ci={qx,Qy) denotes the planar momentum 
operator, and cr= (cti, (72, 173) is the Pauli spin matrix. 
V{r)— UicroS (r— fi) is the potential of two nonmag- 
netic impurities located at ri= (— (i/2,0) and r2— (d/2,0) 
with strength Ui. do is the 2x2 unit matrix. 

The features we discuss are expected to be seen in the 
change of the real-space LDOS owing to the influence of 
magnetic flux which passes through the area enclosed by 
the two scattering paths shown in Fig. [TJ This quan- 
tity can directly reveal the WL or WAL effect in TI via 
AB oscillatory periods in LDOS. The real-space Green's 
function involving the impurities scattering is given by 
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FIG. 1: (Color online) Surface electronic interferometer, 
comprising a spin polarized STM tip at r and two impuri- 
ties at ri and r2 separately. The interference contributions 
in the LDOS is introduced by the electrons traveling along 
clockwise and anticlockwise loops enclosed by the STM tip 
and two impurities. The applied magnetic field afi'ects this 
interference via the AB effect. 



Dyson equation G^Gq+SG, with 

(5G = y" dr"Go (r - r"; w) V (r") G (r", r'; lu) . (2) 

The bare electron Green's function Go for the TI surface 
states could be expressed in terms of a complex amplitude 
multiplied by a "complex" spin rotation: 



Go (r,r';a;) 



-D+R{a,l3,ri) , 



(3) 



which is useful to gain more physical insight into the 
transport between r and r'. Here, D-^-—^/g^+gf, where 
5o/i=sgn(w+)ro/i (|fc|p)T«^o/i (I^Ip) with p=|r-r'|, 
k=kp {l+uj/ Ep), Lo^—Lu+Ep, and Ep=vkp the Fermi 
energy of TI. Jq/i (Vq/i) ^re the first (second) kind 
Bessel functions. R{a, /3,ri) —e^°"^^ e''^°'^e'^^'^^ is a spin 
rotation operator characterized by the three Euler an- 
gles a=-|-f-f with e'^ = {p-x+ip-y)/p, /3= tan"! (|^), 
and r]=—a. Equation ([3]) is exact under a high-energy 
cutoff [33,|3l|. 

Following the perturbation approach, Eq. ([2|) can be 
expanded to any order in the impurity potential Ui . Our 
effort will be concentrated on the scattering processes of 
surface electrons with the both impurities, in which the 
scattering paths enclose loops. Therefore, taking all this 
into account, after a long algebra calculation, we have 

SGl - Go (r - n) I^iGo (ri - rs) T2G0 (ra - r') 

+ Go (r - ra) W2G0 (ra - ri) TiGo (ri - r') , (4) 



where 



W,/2 = 



Tl/2 



c^o - T1/2G0 (ri/2 - r2/i) T2/1G0 (r2/i - ri/2) 



effect from the scattering with both two impurities par- 
ticipated. In the absence of SOI, the two terms in 
Eq. (H}, i.e., the scattering amplitudes corresponding 
to the time-reversal processes, should be equal to each 
other and thereby give rise to the constructive quan- 
tum correction to the LDOS in the so-called WL the- 
ory. In the presence of a strong SOI, however, the 
electrons interference is affected remarkably due to the 
non-collinear multiple scattering trajectories that gener- 
ate nontrivial spin rotations. The spin rotation leads to 
a destructive interference in LDOS by a phase change 
picked up during the clockwise and anticlockwise scat- 
tering processes, which is the origin of the WAL ef- 
fect in systems with strong SOI. Explicitly, for collinear 
scattering paths between ri and r2, there is no net 
spin rotation since Go (ri/2-r2/i) Go (r2/i-ri/2) occtq. 
Whereas for non-collinear multiple scattering tra- 
jectories, such as the loops shown in Fig. 1, 
Go (r-ri/2) Go (ri/2-r2/i) Go (r2/i-r) is not a unit 
matrix which implies net spin rotations during scattering 
processes. As a result, in contrast to the collinear scat- 
tering process, the non-collinear multiple-impurity scat- 
tering on TI surface can induce dramatic modulations in 
the LDOS. 

Applying a magnetic field tends to destroy the destruc- 
tive interference and may bring forth some key signatures 
of WAL effect, such as the $o/2 periodic AB oscillations 
in the spin components of LDOS (similar to $o/2 oscil- 
lations in the magneto-conductance due to WAL). In the 
presence of a low magnetic field, the Green's function can 
be semiclassically approximated as [33| . 



Go (r - r') = e'^ ^('^'^'Go (r - r') , 



(6) 



where A=: (—By, 0,0) represents the vector potential. 
This approximation is exact so long as the magnetic 
length is much greater than the Fermi wave length 
[Ib= ($o/27rB)^/^ >Af]. For magnetic field B=5~10 T, 
the corresponding magnetic length Zb~11-63~8.22 nm, 
while the Fermi wave length Ai?=7.21~5.3 nm for TI with 
i?_F=250'^340 meV Accordingly, the condition in Eq. 
^ for the semiclassical approximation is valid for these 
low magnetic field and low energy ranges. Also, the Zee- 
man splitting is negligibly small (typically of 1.0 meV at 
5=10 T for Bi2Se3 film) compared to the strong SOI, and 
thus is not taken into account in the following discussion. 

The correction of the LDOS due to the magnetic fiux 
is given by 



ANl [r,uj,B) = -i^Tr 

TT 



6Gl(v,^)~5Gl{v,uj) , (7) 



(5) 



where 5G l is calculated from Eq. ^ with Go . For large 
is diagonal with T-matrices Ti/2=j3^j^^;^^|^q:^. Equa- distances (|A:|p>l), Jo/i«±^2/ (tt p) cos [7r/4=F |fc| p] 



tion (U) is a general formula describing the interference and yo/i~±\/2/ (tt |fc| p) sin p=F7r/4]. Combining Eqs. 
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dHZl), we get 

ANl (r, w, B) = AnI (r, oj, B) + AA^^ (r, w, B 



/27r$\ 


- 1 


cos 







(8) 



where Cocfc-^/^y^cos (^) sin (^) sin (^) 
with e'''i/2=(pi/2-x+ipi/2-y)/pi/2 and Pi/2=r-ri/2, 
X(r)=fc(rf+Pi+P2). ^ith 
t{io)~ yiJ^ : is the nonzero element of 



2x)^ (ic^)^-(„<,)^ 

T— matrix. One can notice that [^q [d) —gf (d)] vanishes 
within the above asymptotic representation of Jq/i and 
Yo/i- 

In the present setup, we focus solely on the (spin- 
resolved) LDOS at r={x, y), which is probed by the STM 
tip with the same plane coordinates {x,y). Since the 
STM tip also dually plays as a scatter for the electron's 
loop motion, thus, in the presence of a fixed magnetic 
field B, the AB oscillation varies with the tip position 
along the y direction. Obviously, the AB oscillation sig- 
natures result from cos (^^) — 1=0 in Eq. (|8]), and the 
spacial oscillation period is do—^^ for fixed B and fixed 
impurities configuration, corresponding to a $o period 
in the scale of flux, which gives rise to WL effect in the 
total LDOS. To observe a complete AB oscillation pe- 
riod in ANl (r,w,i?), the magnetic field should satisfy 
the relation 



B2^ - 2^o/yd. 



(9) 



Now let us analyze the spin-up and spin-down compo- 
nents of LDOS. The extraordinary SOI in TI affects the 
surface electron interference due to the non-collinear mul- 
tiple scattering trajectories that generate nontrivial spin 
rotations, which can be easily obtained by the rotation 
operator R{a,/3,ri) in the Green's function Gq. Explic- 
itly, the AB oscillations of spin components of LDOS are 
given by 



AN^^ {r,uj,B) « TC.Im [f(tj)e^(xW-f ) 



(10) 



27r$ 

sm I — — T 9 ] ± sin ( 
<3>n 



where (/)=2L^ 



5[tan-H^^)-tan-i(^)] and 
Cs = Y sin0. As clearly seen from this equation, compar- 
ing to the total LDOS, there occur in the spin-resolved 
LDOS additional strong SOI induced quantum interfer- 
ence signature. The above equation should be reasonable 
because of the destructive interference brought about by 
the spin rotation during clockwise and anticlockwise scat- 
tering processes as well as by the magnetic field. This 
strong spin interference effect deviates the real-space AB 
oscillations from dp period when a SP-STM tip scans on 
the TI surface in the presence of a fixed B, which is our 



concentration in the present paper. We can easily find 
that the spacial AB oscillation period of ANj^ is de- 
termined by the factor ^'1'/-'-= sin(^^=F(/))± sin0. Tak- 
ing AA'^^ as an example (Similar analysis can be done 
on AA^^) to find the AB period, we consider the roots 
of ^^=0 as an equation of y, which can be rewritten 
as Fl~=sin0[cos(^)-l]-cos(/)sin(^)=O. There are 
two cases that result in F'^—O: (i) First, it is easy to get 
that 2/1=^^ (n=0, ±1, • • •) are the roots of ^^^^=0; (h) 



Bd 

Second, expanding F^^ at point y2 — '^^"g^'''^" ; a simple 



equation for y can be obtained M{y 



_ (2n+l)$o 
Bd ' 
(>i+l)*o - 



-P=0, 



Bd 

which gives out other asymptotic roots y= +^ 
{M and P are the expanding coefficients). These roots 
lie near y2 and become better for the limit of small B [see 
following Figs. EJb) and[2ljc)]. Combining with case (i), 
the AB oscillation signals for AA''^ occur at with a 

spacial period of -^ji—^ (i.e., ^ in the scale of flux), the 
half of ANl ■ This half period could be understood as an 
analog of WAL effect in the spin-polarized LDOS since 
strong SOI in TI can result in the two-dimensional WAL 
effect in the magneto-conductance, which can be repre- 
sented as the AB oscillation phenomenon with period of 

$0 

2 • 

Moreover, the interference signature of LDOS decays 

^ /2 

by a factor of [(ipip2] in the asymptotic representa- 
tions, Eqs. ([U and (|10l) . Actually, however, dephasing 
processes have been observed in tr ansp ort investigations 
in Bi2Se3 and Bi2Te3 films [il,|IlH2i as well as in AB- 
effect studies of Bi2Se3 nanowires The phase 

coherence length of Bi2Se3 and Bi2Te3 can be as large 
as hundreds of nanometers, which is tens times of the 
Fermi wave length. The characteristic distance in our 
setup must be much smaller than the phase coherence 
length {d<lcj,), so that we choose d=2Onm^Z0 without 
taking into account the dephasing processes in the fol- 
lowing numerical calculations. 

To supplement the above analytical results, typically, 
we present our numerically calculated data in Fig. [21 
where the upper (lower) panels correspond to B=5 T 
(10 T). The blue horizontal strips are the AB oscillation 
signals in the real-space LDOS, while the oscillatory el- 
lipse features are the interference signals arising from the 
contributions of Nl (r,cj,_B=0). It is obvious that in the 
case of B=5 T (10 T), the interstrip distance is do^Sb 
nm (42 nm) in the total LDOS as shown in Fig. [21(a), 
corresponding to the $o period of AB oscillations. The 
numerical data in Fig. ^a.) are exact while Eq. ((Sj) is 
approximate. 

However, differing from the <I>o periodic AB oscil- 
lations in the total LDOS, Fig. Mh) and W[c) show 
that the interstrip distance in the spin-resolved LDOS 
is corresponding a ^ period. This analogous phe- 
nomenon of WAL with a period of ^ in AA^2^^ i^^ '^j B) 
found from the numerical calculation are accord well 
with our analytical result given by Eq. (jTUj). As ad- 
dressed above, the present analogous WAL phenomenon 
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FIG. 2: (Color online) Simulations of the AB oscillations of 
the electronic LDOS in Bi2Te3(lll) surface with fixed mag- 
netic field B—5 T (upper panels) and B—10 T (lower pan- 
els), (a) the total LDOS patterns; (b) the spin-up component; 
(c) the spin-down component. The parameters are chosen as 
v=2S7 meV-nm, £;f=250 meV (Af=7.21 nm) and d=20 nm. 
The white dots denote the positions of adatoms. The pe- 
riodic horizonal blue strips in patterns are signature of AB 
oscillations in LDOS. 



can be understood via the Berry phase. Explicitly, 
the quantum phase difference between the closed loops 
in opposite directions in our setup corresponds to the 
Berry phase associated with spin rotation by 2tt, which 
is given by Aip=-i Jg"'^'^'' dO {ipq\ ^ \tljq) =±7r. Here, 

\ipq) =^[±.ie^^^'^\l)^ are the eigenstates for the free 
Hamiltonian of TI surface with taiv[9 [q)] . To ex- 
perimentally verify our predicted AB interference strips 
with ^peroid in the spin-resolved LDOS shown in Fig. 
[i;b) and[2i;c), the SP-STM tip and sufRciently strong 
scattering potentials are required, which we believe are 
achievable in current experimental capabilities. So we 
hope the present prediction can be directly observed by 
STM in situ measurement instead of the complicated 
low-temperature transport measurement. 

For comparison, we also calculate the AB oscillations 
on the conventional metal surface with weak but observ- 
able SOL We choose Au(lll) as an example, in which 
the Rashba SOI is ^^40 meV-nm and the Fermi wave 
length is ^3.74 nm [s^. The consequent calculated 
results are shown in Fig. [3]for _B=10 T. From Fig. |3]the 
SOI influence on the electron interference in Au(lll) 
can be summarized as follows: (i) The AB interference 
strips in the spin-resolved LDOS display oscillatory 
behavior along x axis, which is caused by the spin 
rotations induced by non-collinear multiple scattering 
trajectories; (ii) The SOI destroys the elliptic features 
in the LDOS maps even if there is no applied magnetic 
field, see the longitudinal extending blue strips in Fig. |31 
However, the AB oscillation period in the LDOS pattern 
is also <I>o (the interstrip distance is ^42 nm); (iii) 
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FIG. 3: (Color online) Simulated STM maps of the AB oscil- 
lations of the LDOS for two identical impurities 20 nm apart 
on the Au(lll) surface with a SOI constant 7=40 meV-nm 
under B=10 T. (a) the total LDOS pattern, (b) and (c) the 
spin-up and spin-down components, separately. The effec- 
tive electron mass m*=0.26mo and the Fermi wave length 
Af=3.74 nm have been chosen. 



Especially, there is no $o/2 period in the spin-resolved 
LDOS in Fig. [3] because there is zero net Berry phase 
and no topological chirality in the Shockley surface 
state on Au(lll), which is totally different from the 
case of TI surface. Analytically, the unperturbed spatial 
Green's function for the conventional metal surface 
with weak Rashba SOI described by the Hamiltonian 
Hq'^'^ = {hq)'^ /2m*— {'j/h) {crxq)-z is asymptotically ex- 
pressed as Go(r,r',w)«T[/^(p,w)+/i,(p,u;)(<T-p)], 

where r=— j-^tJ'} ^ , p=r r' , and 



_2(fei + fe2)?i2v^' 
P ' VP 



with 



L2 



2m'{Lj+Ep) 



and kso^m*j/h?. After 
a tedious derivation, we find that for the conventional 
metal surface with the same interferometer setup shown 
in Fig. [U the correction in the spin-resolved LDOS due 
to the magnetic flux can be approximated by 



AN 



(r, cv, B) « -2 Im [t (c.) t^/m (Pi) /m {d) f, (P2)] 
X [cos(27r$/$o) - 1] , (11) 



from which the absence of ^ AB oscillatory period be- 
comes obvious. 

In summary, we have performed a semiclassical analy- 
sis of the SP-STM probed AB oscillations in the LDOS 
induced by two impurities on a TI surface as well as on a 
conventional metal surface with spin splitting. We have 
found that the total LDOS in both systems present WL 
phenomenon with an oscillatory period $0 in the AB os- 
cillations. Remarkably, the analogous WAL signified by a 
oscillation period has been found in the spin-resolved 
LDOS in the TI system, while it was absent in the con- 
ventional metal surfaces. This phenomenon, which can 
be observed in the SP-STM experiments, may provide an 
important signature for the existence of the topological 
surface states and provide a useful criterion to distinguish 
the TI surface from other two-dimensional systems. 
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